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ABSTRACT 



In 1997, an extension of general relativity was proposed that predicts the 
dark energy density A to vary linearly with the total number of macroscopic 
black holes in the entire three-volume of the universe. We explore this prediction 
and find that A must be roughly constant after the bulk of the stellar mass black 
holes are in place, so for a redshift z smaller than unity. Conversely, the highest 
black hole formation rate corresponds to the peak in the cosmic star formation 
history, earlier than z = 1. This yields a fast declining A, by a factor of about 
5, from redshift 1 to 3. At even earlier times, before many stars were formed, 
the value of A should be much smaller than its current value. These predicted 
effects are all consistent with current data, and near future observations can 
definitively confirm or disproof the link between the dark energy density and the 
total number of black holes in the entire universe. 



Subject headings: Cosmology: theory — dark energy — black holes 



1. Introduction: Dark Energy and Theories of Gravity 

The expansion of the universe appe ars to be accelerating, as indicated by the magnitude 



redshift relation of type la supernovae (lAmanullah et al. 



Garnavich et al. 1998 



Schmidt et al.l Il998t iRiess et al.l 119981 ) . This phenomenon is usually 



2010; Perlmutter et al. 1999, 1998 



referred to as dark energy and constitutes one of the big mysteries in cosmology. Dark 
energy seems to be fundamentally connected to the dynamics of space-time itself. The 
simplest expression of dark energy is in the form of a positive cosmological constant term on 
the left hand (geometric) side of the Einstein equation. The origin of such a constant is not 
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specified by general relativity (GR), though, leaving one with an incomplete description of 
space-time dynamics. 

The success of Einstein gravity is overwhelming. Hence, any alternative theories should 
take great care in preserving the beautiful features of GR, while pursuing the mystery of 
dark energy. Often, creative and compelling ideas notwithstanding, ad hoc elements or 
unconstrained parameters and fields are part of modified gravit y theories, e.g., vector-tensor 



theory, quintes s ence, M-theory, Einstein-aether theo r y, MOND (IJacobson fc Mattingly||2QQl 



Jacobsonll2010l ; IZlosnik. Ferreira. fc Starkmanl 12007 



Hofava 200 



9; 



Sanders 



2011 



An extension of GR has been proposed in ISpaanj (119971 ). S97 from hereon. This theory 
leads to a unique and testable prediction for dark energy, while fully preserving GR, without 
the introduction of any additional degrees of freedom. The S97 paper is somewhat techni- 
cal (it uses the mathematics of algebraic topology) and we therefore provide an extensive 
summary of its salient physical properties below. 



2. Physical Aspects of the S97 Theory 

The Einstein equation has many symmetries, but it is not invariant under local con- 
formal (scale changing) transformations. As a result, one finds larger metric fluctuations 
when one goes to smaller spatial scales. In fact, at the Planck scale of /p ~ 10 -33 cm, these 
fluctuations in space-time become huge (even singular) and space-time should take the form 
of a quantum foam with Planck mass fmp ~ 1Q~ 5 si black holes (BHs), so-called mini BHs, 



popping out off and into the vacuum (j Wheeler! 1 19571 ). Space-time thus enjoys the presence 



of mini BHs at the Planck scale, i.e., is multiply connected. A fundamental problem with 
the quantum foam is that it fluctuates violently on a Planck time of tp ~ 10~ 43 ' 5 sec and is 
therefore difficult to conceive as a stable structure. It is important to realize here that GR is 
a local geometric theory. Hence, it allows for freedom in the (global) topology of space-time. 

S97 seeks stability for the quantum foam through GR, which is the natural low energy 
limit of any quantum gravity theory. Specifically, the Schwarzschild solution for a macro- 
scopic BH can be used as the low energy limit of the quantum foam in a topological sense. 
That low energy limit that allows for changes in the geometry of space-time but 

one that also assures all local observers to count the same number of event horizons. I.e., 
the connectivity of space-time, so the topology of the event horizon, does not change if one 
makes a BH larger. Furthermore, such a countable quantity like the number of BHs is invari- 
ant under continuous deformations (homeomorphisms) of space-time, which is the essence of 
topology, unlike geometry, which deals with differentiable transformations (shape preserving 
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diffeomorphisms) . 

This kind of a low energy limit uses the fact that all BHs can be described geometrically 
by just their mass, charge and spin, while they all possess a closed surface as well (the 
event horizon ) , irrespective of their local geometry. Furthermore, BHs have a black body 



temperature (lHawkingllf975l ) and small BHs are much warmer than large BHs (the effective 
temperature goes as 1/mass). However, mini BHs in the quantum foam are transient in that 
they are continuously created and destroyed. Hence, if GR is to be a proper low energy 
limit of the quantum foam (with a well defined number of event horizons for all observers), 
the existence of mini BHs must somehow be subject to the presence of macroscopic BHs. 
Afterall, only the latter can exist (much) longer than a local Hubble time under Hawking 
evaporation. 

In S97, the quantum foam stability problem that plaques GR at the Planck scale is solved 
by allowing only macroscopic (long lived) BHs to induce mini BHs in every local Planckian 
volume ~ Zp. One needs to consider distinct Planckian volumes here since quantum physics 
forces space-time to be quantized at the Planck scale if one considers that observers cannot 
measure with a better accuracy than ~ Zp. The S97 argument for the induction of mini BHs 
is simple: Mach's principle, as interpreted by Einstein, dictates that the global distribution 
of matter determines local geometry. Wheeler's local quantum foam, with its non-trivial 
topology in the form of mini BHs, is also a direct expression of GR. Therefore, following 
Einstein's ideas for geometry, matter that is globally in the form of macroscopic BHs must 
determine the Planck scale topology of GR. 

Topology, with its invariance under continuous deformations of space-time, cares only 
about a countable quantity like the number of BHs in the entire three- volume of the universe 
(like the number of handles on a mug). Hence, the number of macroscopic BHs in the 
entire universe, a global quantity, determines the local occurrence of mini BHs. To be more 
quantitative about this, one must compute how many mini BHs Wheeler's quantum foam 
enjoys in every Planckian volume, for a specific number of macroscopic BHs. I.e., one needs 
to assess with what, possibly time dependent, probability 5 these induced mini BHs are 
created in a Planckian volume. The key aspect of 5 is therefore how it depends on the total 
number of BHs in the entire three-volume of the universe. 

Following S97, one first notes that any topologically non-trivial space-time manifold can 
be constructed by adding together three-dimensional building blocks called primes. These 
primes (e.g., handles, three-tori, three-spheres) lend their name from the property that 
they themselves cannot be written as the connected sum of smaller three-dimensional units. 
Furthermore, interactions among those primes can never change their intrinsic properties 
(like the number of closed loops they contain). Finally, following the superposition principle 
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of quantum mechanics, any property of the quantum foam must be carried by the linear sum 
of these primes. The topological extension of Mach's principle above then states that the 
total number of macroscopic BHs in the universe determines the properties of the quantum 
foam. Macroscopic BHs are primes (half BHs do not exist) and have the same intrinsic 
topological properties as mini BHs in the quantum foam. These macroscopic BHs are a 
linear measure of the global space-time topology, as the connected sum of all those BHs, just 
like the mini BHs in the quantum foam are such a linear measure on the Planck scale. Hence, 
5 is linearly proportional to the number of macroscopic BHs in the entire three-volume of 
the universe. 

This, currently tiny, probability S for a single Planckian volume to enjoy the temporary 
presence of mini BHs is computed below, but its crucial property is that it is linear in the 
number of macroscopic BHs (so counting event horizons) and independent of geometry (so 
volume). The latter property forces all local Planckian volumes in the quantum foam to 
respond as one to the presence of macroscopic BHs, irrespective of whether the universe 
is expanding (and thus adding Planckian volumes). I.e., dark energy is purely topological 
in nature and every observer witnesses the same local topological changes in the Planckian 
quantum foam (so space-time has uniform statistics), subject only to how the low energy 
GR limit drives the formation of macroscopic BHs globally. 

One may think that this almost instantaneous response (within a time slice) of the 
dark energy everywhere to BH formation leads to causality violations, but this is not the 
case: In S97, the equations of motion for the topology of the universe are derived. These 
equations contain no spatial information and are first order in the discretized time variable 
t = ntp, for integers n > 1 that distinguish different time slices of thickness tp. I.e., for any 
time t\ labelling the complete three-space of the universe, the solution one Planck time later 
only depends on the global topology at time t±. Also, under a continuous transformation 
one can always deform three-space in such a manner that any two three-space points are 
brought within each other's vicinity. Subsequently, information on the creation of any BH, 
which merely amounts to raising the count of event horizons by one, can thus be brought 
into the vicinity of any local observer as well. Furthermore, this freedom of continuous 
three-space deformation assures that no such information needs to be exchanged faster than 
the speed of light to inform every local Planckian volume, no matter how distant from the 
BH, of the topology changing event. After all, the distance between the two aforementioned 
three-space points can be made as small as a Planck length by the deformation. One should 
realize here that BHs and observers are not moved physically during such continuous de- 
formations. Rather, it is space-time itself that deforms and facilitates the transfer of BH 
number information through topological identifications between space-time points. 
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The latter can best be imagined as loops that provide short-cuts through four-space 
between any two points in three-space. These loops can be constructed by taking a line 
segment (a path) and identifying two points along it to create a loop as an alternative route. 
Doing this for three dimensions, one constructs a three-torus where each surface of the cube 
connects to the opposite one through the fourth dimension (time). Along these short-cuts 
information travels at maximally the speed of light, obeying local Lorentz invariance, but 
across arbitrary three-space distances. Creating these short-cuts through four-space between 
all pairs of Planckian volumes allows topological information to be exchanged causally within 
the same given time slice of thickness tp. The resulting structure is a lattice of three-tori 
(see S97 for further details). Hence, wherever a BH is formed in some time slice its effect on 
the dark energy density in that time slice is transferred purely topologically through these 
identifications, without the need for causality violations and expressed solely by an update 
on the total number of BHs. This testable topological aspect of dark energy augments the 
local geometric approach of GR, and fully preserves the latter, while it assigns countable 
numbers (BH horizons) to entire time slices. Time is relative of course and one may freely 
choose any foliation for the time slices since the number of BHs is a covariant quantity. 

Because topology allows one to identify all space-time points as above, it is the total 
number of macroscopic BHs in the entire universe, so in all of three-space for every time 
slice, that is pertinent to the dark energy density. Therefore, the spatial distribution of 
macroscopic BHs, or the actual volume of our evolving universe, is irrelevant. In all, S97 
finds the dark energy density A to be linearly proportional to the total number of macroscopic 
BHs in the entire three- volume of the universe, N-q^(z), at a given redshift z. In this, redshift 
is a distinctive measure of time in the usual 3 + 1 decomposition for the evolution of the 
three-geometry and three-topology of space-time. I.e., redshift labels different global time 
slices in the geometrical and topological dynamics of three-space. Algebraically, 



and l\ the size of the universe (after inflation, marking the transition to the low energy 
GR limit) when the first BH forms that exists for longer than a local Hubble time. I.e., 
li is frozen in when the quantum foam stabilizes and each induced mini BH contributes 
about nip ~ 10~ 5 g of mass per Planckian volume lp. Equation [T] shows that, for A^bh(^) 
independent of redshift, any local observer concludes that the dark energy density A in g 
cm -3 remains constant as the universe expands. Conversely, the dark energy density scales 
linearly with the total number of macroscopic BHs that are present in the entire three- volume 




■p 



(1) 



where 




(2) 
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of the universe at a given redshift. Induced mini BHs can drive the expansion of the universe, 
so constitute a negative pressure, because their creation as three-dimensional spatial objects 
requires an increase in four-space volume to actually embed them. I.e., every formed event 
horizon is a closed surface to an external observer. 

The absolute value of A is only roughly estimated in S97 due to uncertainty in A^bh today 
and l\. However, the exact value of A is not important here because this paper concentrates 
only on the evolution of the dark energy density with the total number of BHs in the entire 
universe. Of course, the observable universe may be smaller than the entire universe, which 
is relevant for measuring the total number of BHs through observations. If we assume that 
the observable universe is the entire universe, then one can estimate l\ as follows. Today's 
dark energy density is ~ 10~ 29 g cm -3 and the universe presently contains 10 19 macroscopic 
BHs. The latter number follows if one assumes that 1 in 10 3 stars yields a BH, for a Salpeter 
initial mass function (IMF), and that there are 10 11 galaxies with each 10 11 stars. Hence, 
with these numbers ii~2x 10 14 cm. 

The S97 theory provides solutions for the evolution of the very early universe's topology. 
The number of mini BHs increases exponentially during the quantum gravity phase of the 
universe (so before the first macroscopic BH is formed and the quantum foam is stabilized). 
Therefore, mini BHs can drive an inflationary period at these early unstable times, with a 
dark energy density that corresponds to about one mini BH per Planckian volume (5 ~ 1). 



In all, one has the local space-time quantum foam that (I Wheeler! 119571 ) envisioned, but 
globally constrained by the total number of macroscopic BHs in the entire universe. The 
value of Nbh(z) thus controls the temporal variations in the dark energy density, in the spirit 
of Mach's principle, where global properties of the universe determine local ones. 



3. Testing the Connection between Dark Energy and Macroscopic Black Holes 



The universe appears to be approximately flat based on WMAP7 cosmological obser- 
vations (IKomatsu et al.ll201ll ) and A is about 0.734 ± 0.029 in units of the critical density 
p c = 1.88 x 10~ 29 g cm~ 3 . As discussed above, the dark energy density of our universe is 
linearly proportional to the total number of macroscopic BHs in three-space at any redshift. 
Very practically, this means that if the number of BHs in the entire universe at redshift 1 is 
twice higher than at redshift 2, then A in g cm" 3 doubles from z = 2 to z = 1. This paper 
focuses on the putative evolution in A. Therefore, it is only necessary to compute how the 
total number of BHs in the universe changes with redshift, irrespective of what the current 
value of A is. 
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The bulk of the stars, and thus stellar BHs in the universe appears t o be present as early 



as z = 1, when a (1+z) 4 decline ensues in the cosmic star formation rate (IHopkins &: Beacom 



20061 ). This is earlier than current type la supernova (SN) detections probe, hence one 
expects to find an effectively constant A for z < 1. It is relevant here that type la SNe 
occur with a delay of close to 3 Gyr relative to regular star formation. Conversely, given the 
rapidity with which massive stars/BHs are produced around z = 1 — 3, this epoch should 
exhibit a strong decrease in the number of BHs, and thus A, from z = 1 to z = 3. 



These qualitative expectations can be quantified as follows. IHopkins fc Beacoml ( 120061 ). 
their figure 7, derive the comoving type II SN rate density from the star formation history 
of the universe. These authors adopt a type II SN mass ra nge of 8-50 M and the "BG" 



and "SalA" (Salpeter) IMFs of lBaldrv fc Glazebrookl ( 120031 ). We use their z = - 6 data 



for a Salpeter IMF to estimate the evolution of the dark energy density from today's value 
back to higher redshift, in terms of the comoving type II SN rate density dn^/dt (in yr _1 
Mpc -3 ). In this, we assume that stellar mass BHs are produced by type II SNe with some 
constant efficiency e. The comoving type II SN rate density is used here because such an 
expansion-corrected coordinate system allows one to approximate the relative change in the 
total number of BHs that are formed in the universe at a given redshift, as follows. The used 
data, and thus the derived evolution of A(z), pertain to the observable universe only. Hence, 
we assume that the observable universe gives a good statistical representation of the entire 
universe as far as the change in the total number of macroscopic BHs in evolving three-space 
is concerned. We thus seek A(z) normalized to its present value A(0), which corresponds to 
the current number of BHs in the observable universe. 

The ratio A(z)/A(0) can now be expressed analytically, in terms of dns^/dt for observers 
that measure the rate of type II SNe in their comoving volume, as 

AO) _ /* dt(z) e dns^/dt 
W)~ j°dt{z) edn sn /dt 

and is independent of the constant BH formation efficiency e. We perfor m the integrals in 



equat ion [3] for a flat cosmology and the same cosmic time interval as used by IHopkins fc Beacom 



( 20061 ). t(z) = tr — r- We use redshift as a coordinate, but there is no per- 



tinent dependence of our results on coordinate system because of the normalization in[3j Any 
observer, whatever its coordinate system, will find the dimensionless ratio of two numbers 
to be covariant. One may worry about BHs that are formed beyond our horizon. However, 
as explained above, these BHs do contribute to the dark energy density in our comoving 
volume, but without violating causality. At the same time, the right hand side of [3] is em- 
pirical. So, even though we, as observers, have a special perspective on (the history of) 
the universe in our current time slice, we can still derive the change in the total number of 
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BHs, as three-space evolves, if we live in a representative comoving volume. I.e., one that 
enjoys the same time dependence of BH formation as the whole universe. Figure [1] shows 
the evolution of A(z)/A(0) and constitutes a quantitative prediction of the (parameter free) 
S97 model. 

It is found that A(z)/A(0) changes by less than ~ 30% more recently than a redshift of 
unity. Furthermore, a rapid decline in the dark energy density by a factor of 5 occurs from z = 

I — 3. The shading in figure [TJ represents the obser yational uncertainty in de riving the type 

II SN rate from the cosmic star formation history (IHopkins fc Beacomll2006l ). For redshifts 
z > 5 the shading indicates that a significant amount of star/BH formation may occur at very 
early times (but see the discussion section). The value of the equation of state parameter, 




3 4 6 



Redshift (z) 



Fig. 1 . — Redshift evolution of A derived from the type II SN rate data of lHopkins fc Beacom 
(120061 ) . The shading represents the observational uncertainty in deriving the type II SN rate. 



w 



- p (for pressure p and density p) is estimated to be 



1.023 ± 0.09(stat) ± 0.054(sys) 
flAmanullah et al.ll2010t lAstier et al.ll2006f l and -1.06lj££jjj dSullivan et al.lboilh . for z < 1. 
In figure [U we see that A(z)/A(0) is roughly constant and slightly decreasing with increasing 
redshift for z < 1. This is consistent with observational data, because a constant w = —1.1 
corresponds to an allowed change in A(z)/A(0) of about 30%. The latter follows from the 
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conservation equation 



dp/dt 
P 



-3H(l+p/p) = -3H[l+w{z)], 



(4) 



with the Hubble parameter squared H 2 = 8irp/3 (ILinderl 120031 ). If one expands 1 + w as 
1 + w ~ 5w, the relative change in p over a local Hubble time is about (Hp)~ 1 dp/dt = —35w, 
for constant Sw. Indeed, figure [1] shows that A decreases by about 30%, for a constant 
w ~ —11, from z = to z = 1. Furthermore, the sharp drop in A(z) for z > 1 then shows 
that no constant value of w can mimic dark energy behavior f or z ~ — 2 if the total n umber 
of stellar mass BHs controls the vacuum energy density. Also, lAmanullah et al.l (120101 ) point 
out that although — w is close to unity there is room for evolution for z > 1. This is because 
the SN la data o nly weakly c o nstra in dark energy for such early times. Similar constraints 
are provided b y ISeliak et al.l (120051 ) . based on SNe, SDSS and CMB data, as well as by 
Li. Wu. Yu ( 201ll ). based on SN legacy survey data. 



Furthermore, recent WMAP7 work by iKomatsu et al.l ( 120111 ). provides li mits on the 

properties of time d ependent dark energy parameterized by w (a) = wo+w a (l— a) (jChevallier fc Polarski 
200 ll ; ILinderl l2003l ) . where a is the scale factor (a = jtj). In their table 4, these authors 
estimate wq = —0.93 ± 0.12 and w a = —0.38^0^5, for a flat universe. Their 13% uncertainty 
in wq is easily accommodated in our figure (U while the fidu cial range of w n allow s for the 
factor of 5 change in A over z = 1 — 3. Recent results by ISullivan et al.l ( 120111 ) provide 
wo = —0.905 ± 0.196 and w a = — 0.984^|;q97, which are again consistent with our figure 1. 
Indeed, because the bulk of the stellar mass BHs appears not to have formed before z = 1, 
we expect to have w < — 1 for z > 1. Nevertheless, it is also apparent that a modest im- 
provement in future measurement accuracy around z ~ 1 could confirm or rule out the S97 
model. 



4. Discussion 



One can estimate the number of BHs at z ~ 6 from the soft X-ray background (0.5-2 
keV), which is partly produced by accreting black holes at z > 6. This is relevant because 
the H opkins & Beacom (2006) data stop around a redshift of 6. iDijkstra. Haiman. fc Loeb 
(120041 ) estimate an upper limit for the BH mass density of 4 x 10 4 M Q Mpc -3 . This translates 
into 10 14 BHs for the observable universe at z ~ 6 , if the Pop HI IMF is dominate d by stars 



of 10 - 100 M (IBromm. Coppi. fc Larson! 12002 ; lAbel. Brvan. fc Normanl 12002 ). This is 
orders of magnitude smaller than the number of BHs in the observable universe at z — 0, 
Abh ~ 10 19 . Since the bulk of the stellar mass BHs is not assembled at z > 5, stellar origin 
BHs should play no role in the dark energy density at those early times. 
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A top-heavy IMF may aid the formation of BHs. Studies by iHocuk fe Spaana ( 12010 



20 111 ) indicate that star formation within r c = 100(M/10 M Q ) ' pc of a supermassive 



(M > 10 5 M Q ) BH, accreting at 10% of Eddington, leads to a top-heavy IMF with as 
much as 5-10% of all stars formed being more than 8 M Q . This can provide an additional 
boost to the BH production rate at z > 2 when qua sars accrete rapidly, while also being 
consistent with the SN neutrino constraints for z < 1 (IHopkins fe Beacomll2006l ). The total 
amount of gas involved here is likely less than 10% of the total galactic reservoir, given 
that r c is generally smaller than 1 kpc. Still, if obscured infrared and sub-mm bright high- 
redshift galaxies harbor a large (50%) fraction of all star formation in the early universe 
(for z ~ 2 — 3 ), then the peak in the cosmic type II SN rate may be higher and wider 
than expected (jCasey et al.ll201ll ). Subsequently, this would cause a factor of about 2 higher 
A(z)/A(0) than shown by the grey area in figure [T] around z ~ 3. 

Of course, mergers between BHs may suppress the total number of macroscopic BHs, 
and thus the magnitude of A. The merger rate of stellar mass BHs, which form the bulk 
of the BH populat ion in the universe (rather than superm assive ones), in BH-BH binaries, 
is about 500 yr -1 (jO'Leary. O'Shaughnessy. fe Rasioll2007l ). Given that presently roughly 3 
BHs are formed per second in the observable universe, assuming again that about 0.1% of 
all formed stars yield a BH (for a Salpeter IMF), this is unlikely to decrease A significantly 

The present number of BHs in the universe could be declining if numerous primordial 
black holes were created during the quantum gra vity phase of the universe and are evap- 
orating around this time (IPage fe Hawking! 119761 ). Such primordial BHs can survive, and 
dominate by number today, only if they are more massive than Mq ~ 2 x 10 15 g. For M > Mo, 
primordial BHs withstand Hawking evaporation for more than a local Hubble time since the 
latter scales oc M 3 . However, it is likely that the mass function of primordial BHs scales 
more steeply than 1/M. Consequently, the total number of primordial BHs declines with 
cosmic time (S97). This would lead to a strongly diminishing A with decreasing redshift, 
inconsistent with current observational constraints. 

High precision measurements of w, Wq and w a , as cited in the great efforts above, will 
continue to decrease their error bars to the few percent level for z < 1 type la SNe. They can 
thus probe the modest but real decrease in A shown in figure 1 for z — — 1. Furthermore, 
a combination of type la SNe, galaxy clus tering and weak lensing may be quite powerful to 
constrain systematics in this (I Wang]|2012l ). Of course, the bulk of the evolution in A lies at 
z > 1. In this light, it has been shown recently that active galactic nuclei (AGN) can be used 
for high redshift (upto z ~ 4) measurements of dark energy, using reverberation mapping 
(jWatson. Penney. Vestergaard. fe Davis! 1201 ll ). This novel, and very timely, technique is 
particularly well suited to sample the A(z)/A(0) curve for z = 1 — 3, because AGN can 
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be seen to much higher redshift than type la SNe. This method could thus also confirm 
or disproof the linear scaling between the dark energy density and the total number of 
macroscopic BHs. Of course, because the mean density of the universe increases as (1 + z) 3 , 
the relative influence of the dark energy density is less for z > 1 and high accuracy is required. 

The authors thank Andrew Hopkins for sending the type II SNe data derived from the 
star formation history of the universe, and Marianne Vestergaard for discussions on AGN 
reverberation mapping. 
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